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We consider local magnetic moments coupled to conduction electrons with on-site attraction, in 
order to discuss the interplay between pairing and magnetic order. We probe the ground state 
properties of this model on a one-dimensional lattice through pair binding energies and several cor- 
relation functions, calculated by means of density-matrix renormalization group. A phase diagram 
is obtained (for fixed electron density 1/3), from which we infer that coexistence between magnetic 
order and superconductivity is robust, at the expense of a continuous distortion of the magnetic 
arrangement of the local moments, as evidenced by a strong dependence of the characteristic wave 
vector k* with the coupling constants. This allows us to understand some trends of the coexistence, 
such as the influence of the rare earth on k*, as observed experimentally in the borocarbides. 



The advent of high temperature cuprate superconduc- 
tors in the late 1980's singled out the interplay between 
superconductivity and magnetic order. This was followed 
by the discovery of robust coexistence between some de- 
gree of magnetic order and superconductivity in ternary 
and quaternary rare earth compounds jL2iii as well as in 
heavy fermion matter4^ Very recently, a new class of 
FeAs-based superconductors has attracted a lot of at- 
tention due to their (moderately) high critical tempera- 
ture, and new experimental evidence has been gathered 
indicating that superconductivity coexists with a spin- 
density wave state in some members of the fcrropnictidc 
familyi&i In spite of these experimental advances, micro- 
scopic modelling of coexistence between magnetic and 
superconducting orderings is still in its infancy, and con- 
siderable insight should be gained by investigating the 
competition between these two opposing tendencies. 

With this in mind, here we focus on a specific class 
of materials showing this coexistence, namely, the boro- 
carbides, in which the rare earth clement provides lo- 
cal moments (through their /-electrons) , while supercon- 
ductivity arises from phonon-mediated pairing of con- 
duction electrons: singlet superconductors with cither 
antifcrromagnetic or modified ferromagnetic (i.e., spi- 
ral or domain-like) arrangements have been observed 
experimentally^ We assume that pairing of conduction 
electrons can be described by the attractive Hubbard 
models and that they are coupled to local moments 
through a Kondo-like termi^II! The Hamiltonian then 
reads 

H = -* E ( c U„ + H.c.) -tf5> iT n u 

(i.j).cr i 

+ J^S i .<T i , (1) 
i 

where, in standard notation, the sums run over lattice 
sites, with (i, j) denoting nearest neighbor sites, t sets 
the energy scale (we set t = 1 from now on) , U > is the 



attraction strength, and J is taken positive, thus favor- 
ing an antifcrromagnetic coupling between the conduc- 
tion electron spin cr t = J2 a ,/3=± Au^apc^ (cr Q0 denotes 
the Pauli matrix elements) and the localized spin Sjj for 
simplicity, we take S = 1/2. The two competing tenden- 
cies are clear: as J increases, the Kondo-like coupling 
drives the conduction electrons to form singlets with the 
local moments, at the expense of breaking the pairs. 

In order to determine the properties of this model in 
an unbiased way, we consider a one-dimensional lattice 
and resort to the density matrix renormalization group 
(DMRG)ii»i2iI2>I± to obtain the ground-state \tp ) and 
energy, Eo(N e ), where N e is the number of electrons. The 
aspects of competition which will be highlighted here do 
not depend qualitatively on the fact that one-dimensional 
ordering can only be quasi-long-ranged. 

The Hamiltonian (JTJ) is investigated on lattices with N s 
sites (hence 8^= states in the full Hilbert space) and open 
boundary conditions. For the sake of brevity, we restrict 
ourselves to the density n = N e /N s = 1/3, though we 
have also examined other densities, and the results are 
qualitatively the same, apart from half filling. Lattice 
sizes up to 60 sites were used, and truncation errors in 
the DMRG procedure were kept around 10 -5 or smaller. 

The superconducting state is probed with the aid of 
the pair binding energy, 



Er 



2E (N e + 1) - E (N e + 2) - E (N e ) , (2) 



a positive value of which indicates that the state is super- 
conducting, and of the s-wave pairing correlation func- 
tion, 



P s (r) = (ct ;C J T c 1+rT c H _ ri + H.c.). 



(3) 



The magnetic properties of the local moments are probed 
with the real-space correlation functions, 



S""(r) = (Sf S»), n = x,y,z 



(4) 
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Figure 1: (Color online) (a) Pair binding energy [Eq. ©] as a 
function of the inverse lattice size, for (7 = 4, for the values of 
J/U labelling the curves; data are taken for N s = f8, 24, 30, 
36 and 60 sites, (b) Extrapolated (to N s — > oo) pair binding 
energy as a function of J/U for the values of U labelling the 
curves; the error bars are smaller than the data points. 
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Figure 2: (Color online) (a) Spatial dependence of the pairing 
correlation function [Eq. (j3])], in a double- logarithmic scale, 
for a lattice with 60 sites, and U = 6 in all cases but one, for 
which U — 0; (b) average double occupancy (see text) as a 
function of J/U , for different system sizes. 



and their corresponding structure factors, 

S^W^Ee^-^SfSp, n = x,y,z, (5) 



as well as their sums, 

s(k)= E ^( fc )> 



(0) 



where (...) = (-00 1 ■ ■ ■ |V>o)- 

For given values of U and J, we calculate the binding 
energy for N s = 18, 24, 30, 36 and 60 sites, which, when 
plotted as functions of 1/N S , allows for smooth extrapo- 
lations to N s — > oo; the result is shown in Fig.[]Ja). For 
each U, the binding energy decreases as the Kondo cou- 
pling increases, and vanishes at a critical value, (J/U) c , 
signalling the breakdown of superconductivity; see Fig. 
[ljb), which also shows that (J/U) c decreases as U de- 
creases. It is important to check this against the behav- 
ior of other quantities, for consistency. FigureIHa) shows 
the spatial decay of the pairing correlation function, Eq. 
([3]), for U = 6 and, for comparison, one case for (7 = 0. 
(In order to discard edge effects due to open boundaries, 
in our plots of spatially dependent quantities, we place 
the origin at i = 5.) For J = 0, the squares in Fig. [Ha) 
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Figure 3: (Color online) Results for a lattice with 24 sites, 
(7 = 8, and for different values of J/U: (a) Spatial dependence 
of the local-moment zi-correlation function [Eq. ((4])] (b) same 
for the ^-correlation function; (c) wavevector dependence of 
the structure factor; (d) log-log plot of the spatial dependence 
of the pairing correlation function [Eq. Q] 



reproduce the behavior of P s for the attractive Hubbard 
model, which definitely displays quasi-long-range super- 
conducting order ; 15 ! 16 ' 17 by contrast, when U = [down 
triangles in Fig. HJa)] the system is certainly not super- 
conducting, and the correlation function near the chain 
edge is at least seven orders of magnitude smaller than in 
the previous case. In-between these extreme cases, e.g., 
when J/U = 0.3, the pairing correlation function (cir- 
cles) can hardly be distinguished from that for J/U = 0, 
while for J/U = 0.6, the correlations decay as fast as 
when U = 0. This change in spatial decay rate is accom- 
panied by a significant drop in the average site double 
occupancy, d = (l/N s ) X^^^iT^u)' as ^ can be seen from 
Fig. [2jb). These predictions therefore agree with a su- 
perconducting transition taking place at (J/U) c ~ 0.55 
for U = 6, as determined from Fig. [lja). 

In Fig. [3] we correlate the presence of local-moment 
magnetism and superconductivity, in the regime of strong 
attraction (U = 8). For small values of J/U, the 
local-moment correlations are isotropic, i.e., (SfSf +r ) — 
(SfSf +r ), and a spin-density wave (SDW) with period 
2 (peak of S(k) at k = it) is formed; see plots for 
J/U = 0.2 and 0.55 in Figs. G|a)-(c). Fig. E^d) shows 
that in this same regime of J/U, pairing correlations 
are slowly decaying, consistent with a superconducting 
(SC) state. As J/U increases beyond 0.6, one enters 
a ferromagnetic (FM) region characterized by (SfSf +r ) 
displaying ferromagnetic behavior, while (S^S^ +r ) only 
displays short ranged SDW correlations; this picture is 
supported by the magnetic structure factor now showing 
a peak at k = 0, as it can be seen from Fig. (3]Jc). This 
change in magnetic behavior is accompanied by a dras- 
tic change in pairing correlations [Fig. 0(cl)] : they be- 
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Figure 4: (Color online) Same as Fig. [3j but for U = 4. 



come non-superconducting in the FM region. Again, the 
superconducting transition is consistent with the value 
(J/U) c ~ 0.58, extracted from Fig.QIb). 

For smaller values of U, new magnetic phases appear 
between SDW and FM, in the intermediate range of J/U, 
such as ilustratcd in Fig. [4] for U = 4. For J/t/ = 0.05 
and 0.20, magnetic correlations remain isotropic with pe- 
riod 2 (SDW), while the system is still superconducting. 
However, for J/U = 0.4, the system is superconductor, 
magnetic correlations are still isotropic, but now the peak 
has shifted to k = k* ~ 0.67r; therefore, in this case, su- 
perconductivity coexists with an incommensurate spin- 
density-wave state (ICSDW) of the local moments. Fur- 
ther, when S(k*) is plotted as a function of 1/N S (not 
shown), it shows a steady increase with increasing lat- 
tice size, which indicates a true quasi-long-range ordered 
state. For J/U > 0.5, it is useful to consider the cor- 
relation function S(r) = {Si ■ Si +r ) in conjunction with 
the dimer order parameter^ D(i) = (S; • Sj+i), which 
measures the relative orientation of two successive local 
spins. Figure [3] shows that while S(r) oscillates with r, 
D(i) is always positive, indicating that the local moments 
arc in a spiral state; note also that for J/U = 0.5 [Fig. 
Hfc)], S(0) 7^ 0, which indicates that the local moments 
are not in a singlet state, thus no longer isotropic. This 
state is the continuation (to the U > region) of the 
spiral ferromagnetic (SFM) state found for the Kondo 
lattice model (KLM)iA2. Figures 2(c) and (d) respectively 
show that the FM behavior found for U = 8 only sets 
in for J/U > 0.7, and that pairing correlations become 
strongly suppressed above J/U > 0.5. 

We now briefly discuss the magnetic behavior of the 
conduction electrons. Figure [SJa) shows the correlation 
between itinerant and local moments on the same site. 
Below J/U ~ 0.5 (the critical value for superconduc- 
tivity at U = 4), the conduction electrons are not so 
strongly correlated with the local moments, as a result 



Figure 5: Results for a lattice with 36 sites, [7 = 4, and J/U = 
0.55: (a) the local-moment correlation function (Si-Si+ r ), and 
(b) the dimer order parameter D(i) (see text). 
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Figure 6: (Color online) (a) Localized-spin-electron-spin cor- 
relation function as a function of J/U; (b) magnetic structure 
factor for the conduction electrons [Eq. (JSj>] for different values 
of J/U. 



of pair formation. The magnetic structure factor [de- 
fined in a way analogous to that for the local moments, 
Eq. (jGj) ] for the conduction electrons [Fig. [6(b)] shows 
maxima at k = tt, reflecting the fact that unpaired elec- 
trons tend to develop antifcrromagnctic-like correlations; 
nonetheless, the range of spatial decay of these mag- 
netic correlations between the itinerant electrons is al- 
ways much shorter than that of the local moments (data 
not shown). At J/U = 0.5, both cr(0) and S(0) are non- 
zero, so that rotational symmetry breaks down in each 
subsystem. For J/U sufficiently large, most of the elec- 
trons are unpaired, and they follow the magnetic arrange- 
ment of the local moments more easily: the magnitude 
of the local-itinerant correlations increase [Fig. [6(a)] and 
the conduction electrons achieve ferromagnetic-like be- 
havior, as evidenced by the peak of cr(fc) being displaced 
to k = 0, as in Fig. Mjo)- 

Similar analyses have been carried out for other values 
of U, and the results lead to the phase diagram shown in 
Fig. so that the following picture emerges. For weak 
enough Kondo coupling, the paired electrons are super- 
conducting and the unpaired electrons intermediate the 
coupling between the local moments, leading to the non- 
pair-breaking SDW state; the latter is commensurate or 
incommensurate, depending on the strength of U, but 
each subsystem is in a singlet state. It is interesting to 
note that the period-2 SDW is favored in the large-f/ re- 
gion; the same happens for other fillings, which indicates 
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Figure 7: (Color online) Schematic phase diagram for n = 
1/3: SC stands for superconducting, SDW for spin-density 
wave, ICSDW for incommensurate SDW, SFM for spiral fer- 
romagnetic and FM for ferromagnetism (see text). 

that tightly bound pairs can hop more freely (i.e., with- 
out being hindered by the Pauli principle) if the local 
moments are in a period-2 SDW. As J increases, more 
electrons tend to form singlets with the local moments, 
but superconductivity still survives at the expense of an 
adjustment of the SDW wavevector. At some J C (U), su- 
perconductivity is suppressed: above U ~ 2.5, it is ac- 
companied by the breakdown of spin rotational symme- 
try; below U ~ 2.5, the suppression of superconductiv- 
ity takes place within the incommensurate SDW state, 
and rotational symmetry breaks down at a larger Jr(£7), 
within the normal phase. For J strong enough, the sys- 
tem is normal and ferromagnetic. 

We can now make contact of our results with the 
i?Ni2B2C series of compounds, where R is a rare 
earth. Coexistence between superconductivity and some 
SDW is found for R = Tm (ICSDW), Ho and Dy 



(antifcrromagnctism)^ Assuming U tracks the Dcbyc 
temperature, which in turn tracks the inverse ionic ra- 
dius, we can expect U to grow as R varies from Tm to 
Dy. This trend is consistent with Fig. [7] incommensu- 
rate SDW's are favored in the small U region of the phase 
diagram, while commensurate SDW's are favored in the 
large U region. 

In summary, we have considered a model system com- 
posed of local moments coupled through indirect ex- 
change mediated by pairing electrons. The analysis of 
several quantities calculated for a one-dimensional lattice 
through DMRG indicates that superconductivity coex- 
ists with a magnetically ordered local moment state for a 
wide range of parameters. As the coupling J between the 
conduction electron and the local moments increases, a 
superconducting ground state is preserved at the expense 
of a continuous distortion of the magnetic arrangement, 
as evidenced by changes in the characteristic wavevector 
k* . Superconductivity is suppressed for large enough J 
by two distinct routes, depending on the range of U : for 
large [/, by the appearance of a pair-breaking magnetic 
state with broken rotational symmetry or, for small U, 
within an ICSDW state. This model is surely applicable 
to higher dimensions, preserving most of the qualitative 
features discussed here. Indeed, we were allowed to infer 
the qualitative trend found in the quaternary borocar- 
bidc family of superconductors, in which the magnetic 
arrangement depends on the rare earth component. 
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